Spectral Statistics: From Disordered to Chaotic Systems 
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The relation between disordered and chaotic systems is in- 
vestigated. It is obtained by identifying the diffusion operator 
of the disordered systems with the Perron-Frobenius opera- 
tor in the general case. This association enables us to extend 
results obtained in the diffusive regime to general chaotic sys- 
tems. In particular, the two-point level density correlator 
and the structure factor for general chaotic systems are cal- 
culated and characterized. The behavior of the structure fac- 
tor around the Heisenberg time is quantitatively described in 
terms of short periodic orbits. 

The statistical description of the quantum spectra of 
systems which exhibit chaotic dynamics in their classical 
limit, has been conducted mainly along two routes. One 
is to study an ensemble of similar systems, such as disor- 
dered metallic grains, where electrons experience scatter- 
ing by a random potential. In this approach, ensemble 
averaging is a crucial step done at an early stage of the 
calculation. The results of such a calculation apply to an 
individual member of the ensemble, provided the time of 
observation is long enough. The second route is to char- 
acterize the properties of individual systems by means 
of the periodic orbit theory Q. This is possible for a 
system with chaotic dynamics governed by a Hamilto- 
nian that is simple enough, so that the parameters of 
the periodic orbits needed for semiclassical spectral anal- 
ysis can be calculated explicitly. Averaging in this case 
is usually done over an energy interval which consists of 
many energy levels |2j. This approach is very powerful 
in describing the short time behavior of the system, but 
is faced with significant problems when applied to times 
of order or bigger than the Heisenberg time th = /i/A 
or to energies much smaller than the mean level spac- 
ing A 1^. Despite the obvious differences between the 
two approaches it is believed that to a large extent both 
describe the same physics. In this letter we are apply- 
ing results of the first approach in order to extend the 
periodic orbit theory to times close to th- 

The relation between ensembles of metallic grains and 
ensembles of random matrices (RM) [Q is now well un- 
derstood. The supersymmetric non-linear a model 
actually, provides a microscopic justification for the use 
of RM theory in order to describe the universal features 
of these systems. This formalism offers a routine way 
of calculation of a variety of universal correlation func- 
tions for all Dyson pure symmetry RM ensembles and 



for crossovers between them Q . In view of the growing 
interest in applying the super-symmetry approach to the 
investigation of deviations from universality [^|-^, it be- 
comes important to understand the connection between 
the two approaches described above. 

The object that we analyse is the dimensionless two- 
point level density correlator. 



R{s) = A^p{E)p{E + sA)) - 1. 



(1) 



Here p{E) is the density of states at energy E, A is the 
mean level spacing, and (• ■ •) represents ensemble aver- 
aging in the case of disordered systems, or averaging over 
some interval of energy E if an individual chaotic system 
is considered. The universal form of R{s) is especially 
simple in the unitary case. It is the sum of smooth and 
oscillating parts R{s) = (5(s) - [1 -cos(27rs)]/(27r2s2). 
The conventional perturbation theory for disordered met- 
als [0 can provide only the smooth part of R{s) [ pd| . 
The s ^ 1 asymptotics of R{s) in which the oscillatory 
term, non-analytic in 1/s, is retained was recently eval- 
uated in Ref. This result (for s 0) can be still 
presented as a sum. 



Ris) = Rp{s) + Roscis), 



(2) 



of a pcrturbative term Rp{s) and oscillatory one Rosc{s). 
We rewrite the expression for Rp{s) (see Ref. as 



Rp{s) 
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where a = 2 for the unitary ensemble and a = \ for T- 
invariant ensembles. I'(s) is the spectral determinant of 
a classical operator, namely the diffusion operator: 



V{s)=\{A{e,) 



(4) 



Here are eigenvalues (in units of A) of the diffusion 
equation in the grain, and A{e^) is a regularization factor 
which equals for ^ and unity otherwise . Sur- 
prisingly the oscillatory term Rosc{s), which cannot be 
obtained by a perturbative calculation, is also governed 
by the same classical spectral determinant 2?(s). E. g., 
in the unitary case it has the form 



, , cos(27rs)„, , 

Rosc{s) = ^ , ' V{S) 



(5) 



Since I'(s) is purely classical, it is plausible that for any 
chaotic system there exists a classical operator whose 
spectral determinant can be identified with 2'(s). In 
what follows we shall identify this operator for general 
chaotic systems by a semiclassical analysis of relation ^ . 
For the sake of simplicity we shall consider a two dimen- 
sional system which belongs to the unitary ensemble. 

The semiclassical analysis begins with Gutzwiller's 
trace formula which expresses the density of states 
p{E) as a sum over the classical periodic orbits 



PiE) = -r + --- 



det(M;~/)|i/2^ 



(6) 



where p labels a primitive orbit that is characterized by a 
period Tp, action Sp{E), and Maslov phase i/p-, r stands 
for the number of the repetitions of this orbit. Mp is 
the monodromy matrix associated with the linearized dy- 
namics on the Poincare section perpendicular to the or- 
bit. From here on, energy and time will be measured in 
units of A (e = E/A), and h/A respectively. One can 
substitute (||) into (0) and represent R{s) in the form of 
a double sum over the periodic orbits. Rp{s) is given by 
the diagonal part of this sum. Expanding Sp{e + s) up to 
the linear order in s: Sp{e -f- s) ~ Sp{e) + TpS, we obtain 



■p r 



^ |det(M;-/)|- 



(7) 



The traditional way to deal with the above sum is to 
approximate it by an integral: 



E' 



fiTp) 



det{Mp-I)\ 



(8) 



for any sufficiently smooth function f{t). This approx- 
imation, known as the Hannay and Ozorio de Almeida 
(H&OA) sum rule, holds in the limit t oo where long 
periodic orbits which explore the whole energy shell uni- 
formly are considered. In employing it for the calculation 
of Rp{s), the time t should be restricted to the regime 
where it is much larger than the shortest periodic or- 
bits but still smaller than the Heisenberg time th- The 
result associated with it is therefore the universal one 
Rp{s) = -l/27r^s2 which holds as long as s > 1 ||]. Be- 
low we present a more careful treatment of the sum (|^) 
that keeps the non-universal part of Rp{s). 

Let Ap be the eigenvalue (|Ap| > 1) of the monodromy 
matrix Mp. The area preserving property of the latter 
implies that the second eigenvalue of Mp is 1/Ap. Hence, 



|det(M;-/)ri = |Ap| 



fe=o 



(fc + l)Ap 



rk 



and we can rewrite (0) in the form of a triple sum 



Rp{s) 



OC 

^E(fc + i)E;3*pfc' 



(9) 



(10) 



where 



tpk = |Ap|-iA;*^e'^^-. 



(11) 



Using the relation (||) we can determine the spectral de- 
terminant 'D(s) up to a normalization constant: 



V{s) = \UZ{is)\^ 



(12) 



Upon evaluation of the sum over the repetitions in 
Eq. dlQ), the expression for Z{is) takes the form 



l/Z{is)=l[l[e^p[ik + l)cl)itpk)] 



(13) 



p k=0 



where (f>{x) — t^^ ln(l — t)dt. Notice that the normal- 
ization constant TV plays no role in the perturbative part 
of the two-point correlator. We therefore postpone its 
determination. 

Suppose now that all the periodic orbits are very un- 
stable, namely |Ap| ^ 1 for all p-s. In this case tpk 0, 
4'{ipk) ~* ^tpk and Z{z) reduces to the dynamical zeta 
fimction Q, 



= n n ( 1 

p fc=0 
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|Ap|A^ 



(14) 



This function is the spectral determinant associated with 
the Perron-Frobcnius (PF) operator £* (also known as 
Ruelle-Araki or the transfer operator) |l5|. £* is the 
classical evolution operator which propagates phase space 
density for a time t > 0. Its kernel is therefore given by 



^\y,x) = S[y- u{x; t)], 



(15) 



where y and x are phase space vectors representing co- 
ordinates and momenta, and u{x; t) is the point in phase 
space to which a particle that starts its motion at x ar- 
rives after time t. The eigenvalues of the PF operator 
are of the form e""''''*. They are associated with the de- 
caying modes of a disturbance in the density of classical 
particles exhibiting chaotic dynamics, analogous to the 
diffusion modes of disordered system. Yet, the difference 
is that, unlike in the latter case, here 7^-s can appear also 
in complex conjugate pairs 7 = 7'±i7" where 7' > 0. The 
leading eigenvalue of the PF operator, 70 = 0, corresponds 
to the conservation of the number of particles. The dy- 
namical zeta function (|lj) is the spectral determinant 
associated with the eigenvalues 7^. Up to a normaliza- 
tion constant it is given by the product 



l/Z(z)=[]B^(z-7^), 



(16) 



p,fc 



r=l 



where are regularization factors introduced to make 
the product converge. 

Unlike the periodic orbit theory in quantum mechan- 
ics which gives only the leading asymptotics in the limit 



h ^ 0, the periodic orbit expansion of ( [l^ ) is ex- 
act. It is however proper to comment that, in its present 
form, Z{z) cannot be used to determine the eigenvalues 
7^. For this purpose a resummed formula is required. It 
can be obtained by expanding the infinite product over 
the periodic orbits and ordering the various terms in a 
way that leads to maximal cancelation among them. This 
so called cycle expansion exploits the property that 
the dynamics of chaotic systems in phase space is coded 
by a skeleton of few periodic orbits. In particular, the 
long periodic orbits may be approximated by linear com- 
binations of few short ones. 

^From (lej) and (|) it follows that 



1 



27r2 



(17) 



in complete analogy with the result of Ref. ||ll[ for dif- 
fusive systems. The universal part of Rp{s), that was 
obtained using H&OA sum rule, thus corresponds to the 
first term in the sum ( p^ ) (70 = 0). The rest of the sum 
is apparently system-specific. 

We turn now to the determination of the normalization 
constant J\f introduced in (|l^). We shall assume that 
the leading eigenvalue 70 is of unit multiplicity (this is 
the case when the system is ergodic). Comparison of 
Eqs. ([1^) and (||) gives the normalization factor: 



lim zZ(z). 

z~,0 



(18) 



It is customary to express the semiclassical density of 
states as the logarithmic derivative of the Selberg zeta 
function. The latter is defined as the spectral determi- 
nant associated with the semiclassical energy spectrum 
of the system under consideration: 

°° ( piSp{e)-Wp\ 

Ue)^Y[b,ie-e,)=l[l[(l-j^^^], (19) 



p k=0 



|A,|i/2A^ 



where bj are regularization factors, and ej are the semi- 
classical energy levels of the system. The second equality 
above holds for two dimensional systems. One can show 
that the spectral determinant Z{is), satisfies the relation 

Z(is) =exp{(ln[C.(e + 5)]ln[C(e)]}d}. (20) 

where (• ■ ■)d represents an averaging which retains only 
the diagonal elements in the double sum. Since Ap(e) — 
1 — (9/7r9e)31nCs(e -I- iO), the two-point correlator can 
be written as 



Ris) 



-92 



(51nC(e + s)31nG(e)). (21) 



The diagonal approximation gives the perturbative term 



Rp{s) 



-92 



(QlnCs(e + s)31nG(e))d. (22) 



The difference between Eqs. (^ij) and (|2^) can be also 
expressed through the diagonal average. Using Eqs. (|2C) 
and ( p^ ) it is easy to see that Rosds) is given by Eq.^) 
with 



V{s) = exp {25ft (InCs (e + InQ (e))^} . (23) 

It is convenient to present the result in terms of the 
Fourier transform of the two-point level density correla- 
tor, S{t) = / dse^'^'^ R{s) , known as the structure or the 
form factor. RM theory predicts that for the unitary en- 
semble S{t) — min{\T\/2TT,l) (see the light line in Fig. 
1). In the general case 

S{t) = Sp{t) + i [Soscir + 27r) + Sosc{r - 27r)] , (24) 

where Sp and Sosc are respectively associated with the 
perturbative (H) and the non-perturbative (^ parts of 
the two-point correlator. Assuming that the multiplicity 
of all the eigenvalues 7^ is one. 



■7m I • 



(25) 



Again the universal part of S'p(r) associated with H&OA 
sum rule comes from the leading eigenvalue 70 = 0. The 
higher eigenvalues will contribute corrections which are 
in general oscillatory and decrease exponentially. For in- 
stance, the complex pair 7^ ±^7" will contribute the term 
|T-|e-Tikl (2os('^J'r)/7r. The oscillatory part of the struc- 
ture factor can be written as, 



Sosc{t) 



2tt 



E 



^M(!2M)g-7.|r| 

27r7^ 



(26) 



where 2?^(s) is given by 



(27) 



For example, in the case of quasi one dimensional diffu- 
sive system, where the eigenvalues are of the form 7„ = 
Dri^ one can show that 'Dn{iDn'^) = — 4n(— 1)"/ sinh(7rn), 
while for equally spaced eigenvalues 7„ = vn it is 
'Dn{ivn) — l-KJij sinh(7rn). In general it is expected that 
the contribution will come only from the lowest eigenval- 
ues of the PF operator. 

In what follows it will be assumed that the non- 
universal behaviour is dominated by one eigenvalue (or 
possibly a conjugate pair) 71, i.e. 7^^7i for all ii>\. 

In characterizing 5'(t), five domains of the parameter 
T, drawn schematically in Fig 1, are identified: (I) r ~ Tc 
where Tj, is of order of the period of the shortest periodic 
orbit. Here S(t) is composed of several 5-function peaks 
located at the periods of the short orbits and weighted 
according to their instability. (II) Tc < r < 1/7^. De- 
viations from universality associated with (pq) may be 
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FIG. 1. A schematic draw of the structure factor of chaotic 
system belonging to the unitary ensemble. The light line rep- 
resents the universal RM theory result. For the sake of clarity, 
the non-universal features have been exaggerated. 



the general case. This relates the spectral determinant 
associated with the diffusion equation in the grain to the 
dynamical zeta function which can be expressed in terms 
of the classical periodic orbits. We used these relations 
to extend the theory of the structure factor of disordered 
grains to general chaotic systems. It would be interesting 
and important to derive these relations for generic chaotic 
systems. In this respect the recently proposed cr-model- 
like approach for ballistic systems jl^ looks promising. 
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noticeable also in this interval. Their period of oscilla- 
tion 1 / 7" is of the order of Tc ■ (III) l/"f[<T <2Tr - 
the universal perturbative regime where S{t) = t/2tt. 
This is the domain where H&OA sum rule holds. (IV) 
27r— l/Yi<T <2Tr + l/'y'i, the vicinity of the Heisenberg 
time r = 27r. The non-universal features here are in the 
form of exponentially decreasing oscillations very similar 
to those existing in (II). Yet their amplitude and phase 
may be different. In general, the RM singularity at the 
Heisenberg time (the light line in Eig. 1) will be smeared 
out by them. (V) r > 2tt+1/Yj^. Here again the universal 
result S{t) — 1 holds. 

These results can be generalized straightforwardly to 
orthogonal and symplectic chaotic systems: instead of 
Eq. (g) one should use Eq. (5) of Ref. [|| with P{s) — 
s^V{s). 

The behavior of the structure factor in the vicinity of 
the Heisenberg time is a manifestation of a striking prop- 
erty of the periodic orbit sum (H). Namely, that the tail 
of the Gutzwiller's series (the long periodic orbit) en- 
codes its head (short periodic orbit). As a result, S{t) 
in the vicinity of the Heisenberg time is determined by 
the same short periodic orbits as at small r. The argu- 
ment is that the long periodic orbits determine the 
position of the energy levels. Therefore through the long 
range correlation of these levels they encode the infor- 
mation about the short periodic orbits. In fact. Berry 
and Keating resummation method of the periodic 
orbit sum associated with the quantum spectral deter- 
minant of chaotic systems is based on the bootstrapping 
of long periodic orbits with periods near the Heisenberg 
time T ~ 27r to the short ones near Tc- The behavior of 
S{t) near the Heisenberg time reflect this sort of symme- 
try in the sense that it is determined by the short time 
dynamics of the classical system. 

To summarize, we identified the diffusion operator in 
disordered grains with the Perron-Erobenius operator in 
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